A'-slit interference: 
Path integrals, Bohmian trajectories. 
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Path integrals give a possibility to compute in details routes of particles from particle sources 
through slit gratings and further to detectors. The path integral for a particle passing through 
the Gaussian sUt results in the Gaussian wavepacket. The wavepackets prepared on slits 
and superposed together give rise to interference pattern in the near-field zone. It transforms 
to diffraction in the far-field zone represented by divergent principal rays, at that all rays are 
partitioned from each other by (A'^ — 2) subsidiary rays. The Bohmian trajectories in the near-field 
zone of TV-slit gratings show wavy behavior. And they become straight in the far-field zone. The 
trajectories show zigzag behavior on the interference Talbot carpet (ratio of particle wavelength to 
a distance between slits are much smaller than 1 and A'' ^ 1). Interference from the the A'^-slit 
gratings is simulated by scattering monochromatic neutrons (wavelength=0.5 nm). Also we have 
considered simulation of interference fringes arising at scattering on an A''-slit grating of fuUerene 
molecules (according to the real experiment stated in e-print 1001.0468). 

PACS numbers: 03.75.-b, 03.75.Dg, 42.25.Fx, 42. 25. Hz, 45.20.Jj, 47.10. Df, 61.05.fm 



I. INTRODUCTION. 



Classical mechanics operates with point parti- 
cles behavior of which is found unambiguously 
from its variational principles |1]. Initially, we 
begin with the action integral: 



1-1 

S= j L{-^,-^;t)dt 

to 



Here L{'^,'^;t) is a Lagrangian function equal 
to difference of kinetic energy and potential en- 
ergy of the particle. Dynamical variables and 
q are generalized coordinate and velocity of the 
particle. Ones proclaim, that the action S re- 
mains constant along an optimal path of the 
movement particle. It is the least action prin- 
ciple. According to this principle, finding of 
the optimal path adds up to solution of the ex- 
tremum problem 6S = 0. This solution leads to 
Hamilton- Jacobi equation (HJ-equation): 



as 
Jt 



(2) 



Here H{'^,^;t) = (f'^) - L{-^,'^;t) is the 
Hamilton function and = VS* is a particle 



momentum. It is worth mentioning the optical- 
mechanical analogy p] of particular solutions 
of the HJ-equation. It brings to light on deep 
parallels between mechanical trajectories, opti- 
cal rays, and even fiuid streams fall under these 
parallels. 

t 
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FIG. 1: Shift of a surface S = C to a new position S = C + e 
on a value 5 <C 1 [J. Possible trajectories shown by dotted 
blue curves intersect the surfaces perpendicularly. 

Please note, that the gradient of the function 
S, that is, VS, is directed normally to its sur- 
face S = const. Consider two nearby surfaces 
S = C and S = C + e, see Fig. [1} Let us trace 
a normal from an arbitrary point P of the first 
surface up to its intersection with the second sur- 
face at point P'. Hereupon, make another shift 



2 



of the surface on the value of 2e, thereupon on 
3e, and so forth. Until all the space will be filled 
with such secants. Normals drawn from P to 
P', thereupon from P' to P", and so forth, dis- 
close possible trajectories of a mechanical sys- 
tem, since \VS\ = e/6 represents a value of the 
gradient of S. This relation can be expressed in 
the vector form 

t = V^. (3) 

So far as the momentum ^ = m \^ has a direc- 
tion tangent to the trajectory, then the following 
statement is true [1] : trajectory of a moving point 
is perpendicular to the surface S = const. Fig. [T] 
shows possible trajectories by dotted curves in- 
tersecting the surfaces S perpendicularly. 

It is appropriate to mention here the Liouville 
theorem, that adds to the conservation law of 
energy one more a conservation law. Meaning 
of the law is that a trajectory density is con- 
served independently of deformations of the sur- 
face that encloses these trajectories. Mathemat- 
ically, this law is expressed in a form of the con- 
tinuity equation 

1^ + (^Vp) = 0. (4) 

Here p is the density of moving mechanical 
points with the velocity \^ = ~^ /m. 

Thus, two equations, the HJ-equation and the 
continuity equation, give mathematical descrip- 
tion of ensemble of moving points undergoing 
no noise. Draw attention here, that the conti- 
nuity equation depends on solutions of the HJ- 
equation via the term \^ = VS/m. Whereas, 
the HJ-equation does not undergo dependence 
on solutions of the continuity equation. This is 
essential moment at description of moving en- 
semble of the classical objects in contrast to 
quantum ones, as we will see farther. 

In 1933 P. A. M. Dirac drew attention to a 
special role of the action in quantum mechan- 
ics |2] that can exhibit itself in expressions by 
means of a term expliS*//;,}. It is appropriate to 
notice the following observation: the action here 
plays a role of a phase shift. And according to 
the principle of stationary action, the phase shift 
has to be constant along an optimal path of the 



particle. In 1945 Paul Dirac emphasized once 
again, that the classical and quantum mechan- 
ics have many general points of crossing [3]. In 
particular, he had wrote "We can use the formal 
probability to set up a quantum picture rather 
close to the classical picture in which the coor- 
dinates g of a dynamical system have definite 
values at any time. We take a number of times 
ti,t2,t3, ■ ■ ■ following closely one after another 
and set up the formal probability for the g's at 
each of these times lying within specified small 
ranges, this being permissible since the g's at 
any time all commutate. We then get a formal 
probability for the trajectory of the system in 
quantum mechanics lying within certain limits. 
This enables us to speak of some trajectories be- 
ing improbable and others being likely." [3]. 

Dirac's observations had influence to R. Feyn- 
man's searching acceptable language for descrip- 
tion of moving quantum objects, where deci- 
sive role has a term exp{iL6t/h} Idea is 
that the above term proposes mapping a wave 
function from one state to another divided by 
a small time interval St. Feynman's genius in- 
sight has resulted in understanding that the in- 
tegral kernel (propagator) of the time-evolution 
operator can be expressed as a sum over all pos- 
sible paths (not just the classical one) connect- 
ing the points qa and qb with the weight factor 
exp{iS{qa,qb]t)/n} [SIE]: 

K{qh,qa)= A(ix^{\S{qa,qb]T)/h}, (5) 

all paths 

where A is an normalization constant. 

As a result we have now a powerful mathemat- 
ical apparatus, the path integral technique [7]. 
Briefly, according to this idea, there are many 
possible trajectories, that can be traced from a 
source to a detector. But only one trajectory, 
submitting to the principle of stationary action, 
is real. The others cancels each other because of 
interference effects. Such an interpretation is ex- 
tremely productive at generating intuitive imag- 
ination for more perfect understanding quantum 
mechanics. 

The contents of the article is as follows. In 
Sect. [IT] being based on the path integral tech- 
nique we study migration of a particle across a 
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single Gaussian slit. As a result we get a Gaus- 
sian wave packet equivalent to that used in the 



articles P-[TU]. Sect. Ill introduces superposition 
of the wave packets emitting from slits. Here 
we describe also interference from A^-slit grating 
in the near-field zone and diffraction in the far- 



field zone. Sect. [IV] describes the quantum HJ- 
equation and the continuity equation, that lead 
to the guidance equation for finding Bohmian 
trajectories. Simulation of the fullerene molec- 
ular interference stated in P33 is considered in 



Sect. |Vj Sect. |Vl] discusses virtual trajectories, 
virtual sources of spherical waves, and much 
more. Almost all virtual trajectories cancel each 
other. Remaining trajectories are Bohmian ones. 
Sect. 



VII 



concluding section, discusses subtle 
problems of virtual trajectories emerging in vac- 
uum, due to which an optimal paths through 
A^-slit grating are disclosed. 



II. GAUSSIAN SLIT 

Before we will analyze interference on the A^- 
slit grating, let us consider a particle passing 
through a single slit. The problem has been con- 
sidered in detail in [7j. We will study migration 
of the free particle along axis x, see Fig. |2j Its 
Lagrangian is as follows 



X 



m 



(6) 



Here m is mass of the particle and x is its ve- 
locity. By translating a particle's position on a 
small value 6x = (xb — Xa) <^ 1, being performed 
for a small time 5t = (t;, — t^) ^ 1, we find that 
a weight factor of such a translation has the fol- 
lowing form 



, iLSt/h 



exp- 



im[Xb - Xg) 
2h{tb - to) 



(7) 



Pay attention on the following situation: since 
argument of the exponent contains multiplica- 
tion of the Lagrangian L by 5t, putting into 
operation the Lagrangian ^ we obtain result 
{xb — Xa)"^ / {tb — to). Further we will see, that the 
weight factor ([T]) plays an important role. Now, 
by means of such small increments let us trace 



passing the particle from a source through the 
slit and farther. Fig. |2} 



T 
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FIG. 2: Movement of the particle through a slit [7]. 



We suppose, that at the time t = the particle 
leaves the origin of coordinates x = 0. Let we 
know, that after a time T the particle enters to 
the vicinity xq ± 6 of a point Xq, see Fig. |2} The 
question is: what is the probability to disclose 
the particle after a time r at a point xi remote 
from the point xq at a distance Ax = (xi — xq)? 
Let the particle outgoing from the point x = 
at the time t = passes a slit between the points 
XQ — b and x^ + b at the time t = T. Let us com- 
pute the probability of discovering the particle at 
some point Xi after the time r, i.e., at t = T + r. 
Because of existence of an opaque screen, the 
problem can have no solution if we apply laws 
only for moving the free particles along direct 
paths. The particle should pass through the slit 
in order to reach the point Xi. In this connec- 
tion, we partition the problem into two parts. 
Each part relates to movement of the free par- 
ticle. In the first part we consider the particle 
moving from the point x = at the time t = 
to a point x = xq + ^, reaching it at time t = T, 
where |^| < b. The second part deals with the 
particle passing through the point x = xq + ^ 
at the time t = T and moving to the point Xi, 
reaching it at the time t = T + r. A full proba- 
bility amplitude is equal to integral convolution 
of two kernels, each describing movement of the 
free particle: 
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V-lxi, xo) = J K{xi, T + r; xo + e, T)K{xo + T; 0, 0)^^ 

-b 

Here the kernel reads 

K{xb,tb;Xa,ta) = — -exp^ 



m 



im[Xb - Xg) 
2h{tb - to) 



(9) 



It describes a transition amplitude from Xa to Xh 
for a time interval (tb — ta) [Z]. Consequently, the 
integral ^ computes the probabihty amplitude 
of transition from the origin a; = to the point 
Xi through the all possible intermediate points ^ 
situated within the interval (xq — b,xo + b). 

The expression ([s]) is written in accordance 
with a rule of summing amplitudes for successive 
events in time. The first event is the moving par- 
ticle from the origin to the slit. The second event 



is the movement of the particle from the slit to 
the point xi. The slit has a finite width. And 
passing through the slit is conditioned by differ- 
ent alternative possibilities. Therefore, we need 
to integrate along all the slit width in order to 
get a right result. All particles, moving through 
the slit, are free particles and their correspond- 
ing kernels are given by the expression ([9]). By 
substituting this kernel to the integral ([s]) we get 
the following detailed form 



I 



Integration here is executed along the slit width 
a = 2b, i.e., from —b to +b. 

Formally, ranges of the integration can be 
broadened from — oo to -foo. But in this case, 
one should introduce step function G{^) equal 
to unit in the interval [—b, +b] and equal to zero 
outside this interval. Such a form-factor with 
sharp edges can be suitable, if wavelength of the 
particle is much more than interatomic distances 
at slit's edges. If a material from which the 
slit grating is produced shows fuzzy edges of the 
slits, then such a form- factor is incorrect. Sim- 
plest form-factor, simulating such fuzzy edges of 
the slits, is the Gaussian form-factor 

G(0 = eM-e/2b'}. (11) 

Effective width of the curve is conditioned by 
a parameter b. For such a form-factor roughly 
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FIG. 3: Trajectories of particles, passing through the Gaus- 
sian sht ,7,, form a ray with an angle a of the particle beam. 

two thirds of all its area is situated between the 
points —b and +6. If the particles would move 
in classical way, then we should wait, that after 
time r they will be distributed along axis x like- 
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wise as before but with a new distribution center 
Xi be shifted on a value Ax from the point Xq. 
Width bi of the new distribution is also broad- 
ened. The both parameters, Xi and bi, are de- 
termined by expressions 



Xi = Xo 1 



T 



bi = b 1 



T 



(12) 



as shown in Fig. |3j One can see, all classical tra- 
jectories, passing through the Gaussian slit, form 



a ray within an angle of beam a. Observe, how- 
ever, that quantum particles, in contrast to the 
classical ones, at scattering on the slit manifest 
wavelike nature. Relation between matter waves 
and their wavelengths had been established by 
de Broglie. Observe, that the wavelike nature 
underlies phase shifts of the moving particles in 
an observation point. 



Taking into account the above stated, Eq. (10) 
can be rewritten in the following form 



ij{xi,xo) 



T 



Me 



By substituting G{^) from Eq. (11) to this expression and integrating it we obtain 

-1/2 



^{Xl,Xo) =^ 



m 



1 1 



T 



mb'^ 



(hnff {xi - xpf x|\ _ (-(xi - xo)/t + xp/Tf 
^^''^\2h[[ T Tj {l/T + l/T + ih/mb^) 

Here Ax has been replaced by its explicit expression Ax = (xi — Xq). 



(13) 



(14) 



A. Solitary Gaussian slit 



Rewrite Eq. (|14|) by executing insignificant transformations 



T(l + r/T) 1 



m62(l + r/T) 



^1/2 



X exp-! 



im / (xi — xo ) 



2n 



+ 



Xr 



T 



(xi - Xo)/t + Xp/Tfr 



T (1 + r/T)(l + i;ir/mfe2(l + r/T))' 



(15) 



The transformation, resulting to Eq. (15), affects is a plane wave. In accordance with the limit 



only arrangement of the time parameters r and 



T. Now, according to Eq. (12) we could intro- 
duce bi = 6(1 -|- t/T). But in the limit T — )■ oo 
we have equality bi = b. So, we can remove ev- 



erywhere in Eq (15) the term (1 + r/T). Due to 
the limit T — )■ oo we have put source of the par- 
ticles to infinity. According to this trick, we can 
now suppose that the particles incident to the 
slit grating are described by wave function that 



T — i- oo, Eq. (15) transforms to the following 
form 



i!{xi,Xo) 



ml , hr 
1 + i- 



-1/2 



X exp 



27rihT \ mb"^ J 
( (xi-xo)V2&'' 



1 + ihr/mb'^ 



(16) 



Let us clarify correspondence of this formula 
with the Gaussian wave packet, that is described 
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in the articles [8HT0]. This correspondence is 
reached by introducing an effective slit's half- 
width 



a 



y/2 2^2' 



(17) 



where double b, a = 26, is adopted metric 
width of the slit. Now, let us also define the 
complex time-dependent spreading [S] 



a 



. Tit 

[ . 

2ma 



This spreading consists of real and imaginary 
parts. The imaginary part, Axi = hT/2ma, 



emergent in Eq (14), represents uncertainty of 
the particle position on the axis x. What is 
physical sense of this part? First, we rewrite the 
imaginary part in the following form mAxi/r = 
h/2a. Observe that Axi/r is an uncertainty, 
6vx, of a drift velocity of the particle along the 
Consequently, 5px = m5vx is the uncer- 
tainty of its momentum along the same axis. On 
the other hand, let the effective slit width 2(j be 
uncertainty of the particle coordinate, 5x = 2a 
at passing across the slit. As a result, we find 
SpxSx = h. One can see, it is the Heisenberg 
uncertainty principle \7]. 



By replacing relevant expressions in Eq. (16) 



by the parameters ( 17) and ( 18 ) we get the Gaus- 



sian wave packet in the following form 



ip{xi,Xo) = \l 



m 



. pi37r/4 



X exp 



27r^T Y ar 

(xi - Xq)^ 



Aaar 



(19) 



Here we have rewritten the square root of imag- 
inary unit as a factor exp{i37r/4} = 
Pay attention to Eq. (19), the time parame- 



ter T did not disappear fully, but it remained 
in the denominator. In fact, it means, that by 
removing the particle source to infinity, its inten- 
sity, ip'^ip, in the vicinity of the slit grating be- 
come weaker in times. We will suppose, that 
a particle beam from the source is monochro- 
matic, i.e., velocities of all particles are equal 
to the same quantity Vz- Given path length Z 
from the source to the slit grating, the velocity is 



Vz = Z/T. And the momentum of the particles 
is Pz = rnVzi where m is its mass. 

De Broglie relations set momentum pz and en- 
ergy E = p^/2m of the particle in correspon- 
dence with its wave characteristics such as a 
wave vector kz = Pz/T. and a frequency u = E/h. 
In that plane wave mentioned earlier 

can be written down as exp{iut — ikzz}. Now, 
the Gaussian wave packet (19) can be supple- 



mented by this plane wave, that describes, in 
the paraxial approximation [T2], spreading of in- 
cident wave along the axis z: 



ij{xi,Xo,z) 



m 



2nhT 



(xi - Xof 



xexp{i{ut- kzZ + 3Tr/A)}. (20) 

Since in the paraxial approximation r = z/vz, 
then argument of the wave function contains z 
instead of r. In this case, interference experi- 

, 




FIG. 4: Interference experiment in cylindrical geometry. Slit 
grating with n = 0, 1, • • • , — 1 slits is situated in a plane 
{x,y). Radiation goes along axis z. 

ment can be considered on an A^-slit grating sit- 
uated in a plane ( shown in Fig. [4} 
One can see, by comparing equation (20) with 



such a formula presented in pi)J, that with ex- 
ception of some details these formulas are equiv- 
alent. For example, presence of the constant 



phase 37r/4 in Eq. (20) is not critical for further 
observation of interference effects. Essential dif- 
ferences are, however, observed in normalization 
factors. Namely, the term 



A 



ma 
2nhT 



(21) 
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differs essentially from analogous term in the for- 
mula presented in [TU]. Here we see the time 
parameter T ^ 1 is found in the denominator. 
As was mentioned earlier, this parameter has a 
deep sense: as you move farther the source of the 
particles, its intensity on the grating becomes 
weaker. For exception of this effect, the term 
(21) does not exert influence on the interference 
pattern. Therefore, for short we will write this 
normalization factor as A. 

III. ORDERED SEQUENCE OF THE SLITS 

Let a screen, on which a monochromatic beam 
of the particle scatters, has N slits (n = 
0, 1, 2, ■ ■ ■ , — 1) located at equal distance from 
each other, as shown in Fig. |4j Here we have an 
origin of coordinates placed in the center of the 
slit grating. In this frame of reference, n-th slit 
has a position xq = {n — {N — l)/2)d, where d is 
a distance between slits. Further it is measured 
in units multiple to the wavelength A. 

We need now to compute contributions of all 
paths passing from the source through all slits 
in the screen and farther to a point of observa- 
tion (x, z). Per se, we have to execute superposi- 
tion in the observation point of all Gaussian wave 
packets (20) from all slits n = 0, 1, 2, 



N-1. 



Such a superposition reads 



^^^p{x,(n-^yzy (22) 

and probability density in the vicinity of the ob- 
servation point (x, z) is 



p{x,z) = (^(x, z)|^(a;, z)). 



(23) 



Hereinafter, the superscript 1 at x is dropped. 

Before we will take up interference effects in 
the near-field region and behavior of Bohmian 
trajectories here, let us consider an asymptotic 



limit of the formula (23) in the far-field region. 



With this aim in mind, we preliminarily replace 
the term {n - {N - l)/2)d in Eq. by kd, 
where k runs (-(A^-l)/2, ■ ■ ■ , (iV-l)/2). Next, 
at summation we will neglect contribution of 



coefficients at k'^d'^ emergent at decomposition 
(x — kd)"^ = x^ — 2kxd + (kd)^. The point is that 
the terms at k'^d'^ lead messed phases on infinity, 
due to which summation of that exponents gives 
zero contribution. Other sums with coefficients 
at x^ and 2kxd can be easily computed. 

As a result, intensity of the particle beam in 
the far-field region is as follows 



sm 



I{x,z) = Io{x,z)- 



sm 




(24) 



Here terms ({x, z) and /o(x, z) read 
xd - 

C{x,z) = 



47r(7^ 
2a2 ' 



(25) 



Inix, z) 



A^ 



X 



exp< 



2al 



The parameter A is determined by Eq. (21), and 
ay has a form 



\ 



1 



' z\ 
Ana"' 



(26) 



It turns, this parameter is equivalent to the 
instantaneous Gaussian width presented in j8]. 
Emergence of the wavelength A, in the above 



stated formulas, (25) and (26), is due to the de 



Broglie relation binding the momentum with 
the wave vector kz of the pilot wave, namely, 
hkz = h/X = mVz = mz/r. 

Let us now compare intensities of radiation 
from A^ slits in the far-field region computed by 



Eq. (23) and according to the approximated for- 



mula (24). Fig. 5 shows, in gray palette, a dis- 



tribution of the probability density (23) formed 
by radiation from N = 7 slits grating in the far- 
field region. Thermal neutrons are adopted here 
as incident particles on the A^-slit grating. The 
wavelength is A = 0.5 nm. Large red arrows in 
this figure point out directions of radiation of 
principal maxima having a big intensity. Short 
blue arrows, in turn, point out directions of sub- 
sidiary maxima having a low intensity. All near- 
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est principal maxima are partitioned from each 
other by iV — 2 = 5 subsidiary maxima. 




I2Z-J l6Zj 20Zj 

FIG. 5: DifTraction in the far-field zone by scattering thermal 
neutrons (A = 0.5 nm) on A'' = 7 slits grating. Directions of 
principal and subsidiary maxima are pointed out by red and 
blue arrows, respectively. Distance between slits d = lOA, 
width of slits a — 2A, effective width a ~ 0.534a. 




A', m 



FIG. 6: Diffraction in the far-field zone from N = 7 slits 
grating, distance z = 1 m. Circles relate to the probability 



density calculated by Eqs. (22h(23l. Intensity (24 1 is drawn 
by solid curve. Dotted curve draws envelope Io{x, z) ■ . 

Fig. [6] shows diffraction from N = 7 shts grat- 
ing in the far-field zone, distance z = 10^2^ = 
1 m. It is seen, that the principal maxima are 
partitioned from each other by iV — 2 = 5 sub- 
sidiary maxima. The curve depicted by red cir- 
cles relates to the probability density, that is 



calculated by Eqs. (22)- (23). Intensity (24) is Here 



drawn by solid curve. The both curves drawn 
in the figure give a good agreement. The dotted 
curve represents envelope Io{x, z) ■ N^. 

IV. MADELUNG-BOHMIAN INSIGHT 

Majority of known in physics functions are 
representable by linear combination of orthog- 
onal functions, that are solutions of Schrodinger 



wave equation |T3]. Superposition of the Gaus- 
sian wave packets from the slit grating can be 
represented by linear combination of orthogonal 
functions also. So, solutions of the Schrodinger 
wave equation underlie this superposition as 
well. Madelung and next Bohm had shown [HI - 
[T6] . that the Schrodinger equation can be re- 
duced to two coupled equations, the both for 
real-valued functions. With the purpose of the 
following analysis of the interference effect we 
remind these equations. 

Almost at the same time with publication of 
the epochal article of Schrodinger [TB] , Madelung 
published his article entitled "Quantum theory 
in hydrodynamic form" [Tl]. In this article 
Madelung found a deep analogy of behavior of 
the quantum system described by a wave func- 
tion with the behavior of ideal liquid. The first 
step is that, the complex- valued wave function 
can be represented by sum of the functions, 
each is real- valued function. One of this real- 
valued function reduces conditionally to a rule 
for finding velocity of a liquid outflow. The other 
function gives conservation of density of such 
a quantum liquid. Equations, describing these 
two manifestations of the quantum liquid, are 
well known the Hamilton- Jacobi equation and 
the continuity equation mentioned in Introduc- 
tion. They are Eqs. ^ and (|4]). A main dif- 
ference exists, however. The difference is that 
the HJ-equation has an extra term, so-called the 
quantum potential [TSHIT] . 

These renewed equations read 



dS 

Jt 
dp 

at 



= H{q,p;t) + Q{q,p;t), (27) 
+ (vVp) =0. (28) 



H{q, p; t) = + f/(q) (29) 

2m 



is Hamiltonian of the quantum system, where 
f/(q) is a potential energy. The term Q(q, p;t) 
is the quantum potential 



<5(q,p;i) 



2m 



2p 




(30) 
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One can see, that the both equations, (27) 



and (28), are hnked together. Namely, Eq. (27) 



HJ-equation, effects to the continuity equa- 
tion (28) via the velocity v = VS/m. In 
turn, the continuity equation effects to the HJ- 
equation by means of the quantum potential, 
since the latter depends essentially on the am- 
plitude density p. 



Equations (27)-(28) give a rule for computing 
trajectories, Bohmian trajectories. They dis- 
close possible paths of the particles from the 
source to the detector. The rule is given by the 
guidance equation that has a form 

■ (31) 



X 



m 



m 



According to this equation, position (x, z) of the 
particle in a space reaching behind the slit grat- 
ing is given by 



x{t) 







(ir. 



(32) 



_ z(t) = Zq + V-J. 

One could think that the Bohmian trajecto- 
ries are physical artifacts, since they enter into 
a rough contradiction with the Heisenberg un- 
certainty principle, because of prediction in each 
time moment of exact values of coordinates and 
velocities of the particle. However, one should 
remember, that the Heisenberg uncertainty prin- 
ciple relates to procedure of measurements of 
complementary variables. It means only, that 
we cannot exactly measure the both variables, 
the velocity and the coordinate of the particle 
moving along the Bohmian trajectory. 



It should be noted, the computations (32) are 
executed in the paraxial approximation, there- 
fore the coordinate z is determined by a simple 
shift. Fig. [7] shows diffraction of monochromatic 
thermal neutrons (A = 0.5 nm) from N = 7 slits 
grating (distance between the slits d = lOA and 
effective slit width a ~ 0.354a). The diffraction 
pattern shown in this figure covers the near-field 
zone, z < Az^ = 400 nm, see Fig. |5} 

Here and everywhere the Talbot length 

z, = 2- (33) 




FIG. 7: Diffraction in a zone z < 42^ from 7 slits grating 
of thermal neutrons (A — 0.5 nm, distance between slits d = 
lOA, effective slit width a ~ 0.354o). Wavy curves drawn 
by alternating violet and blue colors in the upper part are 
Bohmian trajectories divergent from the slits to infinity. 



is adopted as an important scale parameter that 
divides zones of interference patterns. This 
length starts from Henry Fox Talbot who discov- 
ered in 1836 [18] a beautiful interference pattern, 
Fig|8} that carries his name. 

Fig. [7] shows a probability density distribution 

executed in the gray palette ranging from white 
to black, from zero density to maximal, respec- 
tively. It is clearly seen, that replicas of the grat- 
ing begin to disintegrate already in the nearest 
vicinity of the slits, i.e., at z < 3/2z^. One can 
see, that ai z > 2z^ the replicas disappear fully 
and transition to the far-field zone takes place. 
Namely, transformation of the probability den- 
sity distribution to a characteristic radial radia- 
tion is observed, see Fig. |5] 

Wavy curves shown in the upper part of Fig. [7] 
by alternating violet and blue colors are the 
Bohmian trajectories calculated according to for- 
mulas presented above. Alternation of colors 
permits to discern that the Bohmian trajecto- 
ries do not intersect as the time t, i.e., z = v^t, 
goes on. One can see, all Bohmian trajectories 
diverge as they move away from the slits. Draw 
attention, that the particles, passing through the 
slits, perform undulatory motions similarly to 
colored streams of water leaking across break- 
waters. Such a behavior of quantum objects dis- 
closes analogy with hydrodynamic flows of liq- 
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uids, as it has been described in [HI [19 



1600 3200 ^I'k 4m) 6400 8000 




FIG. 8: Optical Talbot effect for monochromatic light, shown 
as a "Talbot Carpet". The figure has been captured from 
URL http : //en. wikipedia. org/wiki/Talbot_ef f ect 




FIG. 9: Talbot-carpet near slit sources of slit grating, A'' = 
512. Wavelength A = 0.5 nm, distance between slits d = 50A, 
effective width of slits a « 0.354a. In the upper part Bohmian 
trajectories are drawn by alternating violet and blue colors. 

Since we have mentioned above Henry Fox Tal- 
bot and his beautiful pattern, Fig. [8| it would be 
appropriate to consider behavior of the Bohmian 
trajectories on the Talbot carpet [10]. Their be- 
havior on the interference Talbot carpets is im- 
pressive. Fig. [9] shows the Talbot-carpet emer- 
gent by simulation of scattering monochromatic 
thermal neutrons (A = 0.5 nm) on the slit grat- 
ing containing = 512 slits. Distance between 
slits, in this case, is d = 50A = 25 nm, so 
X/d = 0.02 < 1. The Talbot length of that car- 
pet is Zrj, = 2500 nm. And the ratio of this length 
to the distance between slits is z^/d = 100. So, 
such Talbot carpet can be observed in a long 



narrow neutron guide section. 

In Fig. [9]the trajectories, that we see in the up- 
per part of the figure, demonstrate clear zigzag 
behavior. They have tendency to pass through 
caustics (dark patches in this figure) and avoid 
lacunae (white lens-like domains). If a trajectory 
traverses across a light-colored region, it crosses 
this region along a shortest route. 






-4- 




FIG. 10: Interference pattern in the near-field region from 
the grating A*' = 64; A = 0.5 nm, d = 500 nm, z,j. = 1 mm. 
Arrows a and b point out to bohmian trajectories. 
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FIG. 11: Interference pattern in the near-field region from 
the grating A'^ = 64; A = 0.5 nm, d = 500 nm, = 1 mm. 
Arrows a and b point out to bohmian trajectories. 

In conclusion of this section let us consider a 
hard case of the ratio A/rf, namely d = lO^A and 



the ratio is X/d = 0.001. Figs. 10 and 11 show an 



interference pattern in the near-field region from 
the grating, having = 64 slits, and for the case 
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A = 0.5 nm and d = 500 nm. The Talbot length 
in this case is a very large, = 10^ nm = 1 mm. 



We see in Fig. [TOjthe Talbot carpets, which lose 
their forms as z increases. In the vicinity of the 
grating the Talbot carpets demonstrate a hard 
fractality - alternation of the caustics and lacu- 
nae is so dense, that they merge into a gray mass. 
In the limit X/d — )■ 0, in particular, the wave 
function is continuous everywhere, but differen- 
tiable nowhere [20]. If the wave function con- 
tains fractal features, its nondifferentiability for- 
bids a direct calculation of the trajectories. Tun- 
neling processes taking place in those cases mani- 
fest themselves as hopping events. Bohmian tra- 
jectories depicted in this figure show nearby the 
slits such a Brownian-like movement. As a parti- 
cle moves away from the grating, its behavior be- 
comes more predictable - it prefer to pass across 
the caustics and make the round of the lacunae, 
see trajectories pointed out by arrows a and b. 



Fig. 11 shows the same interference pattern, as 
shown in Fig. 
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but here it is extended up to 
z ^ 20z^. We can see, that the Talbot carpets 
exist only within a narrow strip by width ranging 
from z = to about z = iz^. Next, they disap- 
pear, and instead of them a hexagonal packing 
of the interference pattern emerges. It occupies 



a zone from about z 



9z^ to z 



nal stage is transition to the diffraction in the 
far-field zone, z > 16z^. Bohmian trajectories 
drawn in the upper part of the figure behave 
by the Brownian-like manner in the vicinity of 
the slit grating and begin to diverge from each 
other as z increases. Two divergent trajecto- 
ries pointed out by arrows a and b demonstrate 
transition both across the narrow strip occupied 
by the Talbot-carpets and through the hexagon 
zone. In the both cases we see different behavior 
of particles, passing through these zones. In the 
first case, sharp zigzag behavior takes place. In 
the second case, we observe wavy behavior. As 
the particle passes into the far-field region the 
wavy trajectory become more smooth, until it 
turns into a straight line asymptotically. 

All three Berry's conditions for emergence of 
the Talbot carpets shown above are satisfied. 
These conditions read [121 1211423] : (a) paraxial 
beam; (b) arbitrary large number of slits; (c) ar- 



bitrary small ratio X/d. In the limits — ?■ oo 
and X/d ^ the Talbot carpet tends to fractal 
interference carpet. 



SIMULATION OF THE FULLERENE 
MOLECULAR INTERFERENCE 



Atom and molecular interferometry [24] has 
an important significance in modern lithography 
technology. In this key, fullerene molecular in- 
terference experiment was recently presented in 
arXiv.org [TT]. Heavy fullerene molecules Cgo 
by passing through a slit grating disclose inter- 
ference fringes in the near-field zone. Fig. 12 



White dots against blue background visible in 



15Zr^. Fi- o 
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FIG. 12: Fullerene molecular interference. The particle- 
nature (a) and the quantum wave-nature (b) of the surface 
deposited fullerene molecules in one and the same image [llj 



Fig. |12[a) are images of fullerene molecules de- 
posited on a silicon plate. All dots are seen to 
form vertical strips. Summing over each narrow 
vertical band number of the dots we disclose in- 



terference fringe shown in Fig. 12 b). 

Experimental setup and parameters of the ex- 
periment are described in the article in detail. At 
simulating the fullerene molecular interference 
fringes we will use parameters given in this arti- 
cle. Mass of the fullerene molecule Ceo is equal 
to multiplication of the nucleon mass (6 protons 
and 6 neutrons) to 60 nuclei. That is, the mass 
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1.204 X 10-2^ kg. The following pa- 
rameters are presented in [llj: (a) for Cgo with 



IS m„ ^ 

^60 



a mean velocity f , 



111 m/s, a de Broglie 



wavelength X^b is about 5 pm. Observe, that 
it is much smaller then radius of the fullerene, 
RcgQ, that is about 370 pm! The SiN^^ grating 
has a highly accurate period of d = 257.40 nm 
with open slit windows as small as 150 nm for 
the second grating. We take, roughly, d = 250 
nm and a = 150 nm. 

Observe, that diffraction at each of the indi- 
vidual slits within the first grating expands the 
molecular coherence function, due to this it cov- 
ers several slits on the second grating [11]. We 
assume in this key, that the fullerene molecu- 
lar beam has been prepared paraxial coher- 
ent beam [231 - I26] . Simulation of such a fullerene 

k P(x,z) 



m = I.204E-24 kg 


slits = 9 




A = 5£-;2 m 


X-scale — 


1:20 


d = 2.5E-7 m 


Z-scalc — 


1:1250000 


slit width = 1.5E-7m 


P-scale = 


1:2 


a = 2.65E-S m 








FIG. 13: Probability density distribution of the fullerene 
interference simulation. The Talbot length is Zrj. = 25 mm. 

molecular beam incident on the second 9-slits 
grating is shown in Fig. 



13 It is the proba- 



bility density z) of detecting of a fullerene 
molecule within the vicinity of the point {x,z). 
Parameters of the simulation are written in a ta- 



ble inserted in the figure. The Talbot length (33) 
in this case is equal to 2^ = 25 mm. 

At half of the Talbot length, L = Juff- 
mann et al. [26] have observed the fullerene 
molecular interference shown in Fig 



interference fringes, computed at the same 
length, are shown in Fig. [l4Fa). At compar- 



ing Figs. 12 b) and 14 a) we can see their qual- 
itative resemblance accurate to experimental er- 
rors. Collisions between the fullerene molecules 

0.15- 




0.05 



-1500 -1000 -500 500 (ujjj-) 1500 

FIG. 14: Interference fringes observed at some cross-sections 

X, z) shown in the previous figure: 
L = Zt, = 25 mm. 



of the probability density 
(a) L = 2^/2 = 12.5 mm; 



(b) 



occur induced perhaps by dispersion of the ve- 
locities up to ±5 m/s [11]. They distort the in- 
terference pattern, i.e., the pattern get blurred. 
Because of this, the experimental interference 
fringes have no zero level partitions. 

It would be interesting to check the next 
fringes positioned at the cross-section L = z^^. 
They should have maxima shifted on half-period, 
as shown in Fig. [l4Fb). 




FIG. 15: Diffraction in a zone z < 10z,j. from 9 slits grating 
of fullerene molecules beam (A = 5 pm, distance between 
slits d = 5 X 10* A). Wavy curves drawn by alternating violet 
and blue colors in the upper part are Bohmian trajectories 
divergent from the slits to infinity. 

The particle-nature of the fullerene molecules, 
passing through the slit grating, can be evalu- 
ated via the Bohmian trajectories that are cal- 



12 The culated according to the guidance equation (31) 
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and Eqs. (32), where = h/Xm^ 



» 110 m/s. 

Fig. 15 shows contour plot of the probabihty den- 
sity distribution p{x, z) projected on the plane 
(x, z). The quantum wave-nature of the fuUerene 
molecular flow is drawn by color gray. The 
particle-nature presented by the Bohmian tra- 
jectories is shown in the upper part of the fig- 
ure. Predominantly, they are drawn by color 
violet. Some trajectories colored in blue are 
drawn against this background. Due to such a 
color representation flows of fullerene molecules 
through the near-field region and further are eas- 
ily visualized. Fullerene molecular flows are seen 
to have undulatory motions. 

Obviously, there are no possibilities to observe 
real movements of fullerene molecules through 
the slit grating and further [23]. Any perturba- 
tion of the molecule destroys its movement and, 
consequently, it destroys the interference pat- 
tern. In best case we can observe the fullerene 
molecular flows at crossing some kind of detect- 
ing plates (fixed at z^/2, z^, 3z^/2, ■ ■ -, for in- 
stance), as shown in Fig. [I2|a). 



VI. VARIATIONAL CALCULATIONS 

What could force the particle to carry out such 
undulatory and zigzag behaviors, as shown in 
the figures above? Possible answer can be as fol- 
lows: a set of ordered slits in the screen poses 
itself as a quantum object that provides a po- 
larization of the vacuum in the near-field region. 
The polarization, in turn, induces formation of a 
virtual particle escort around of a flying real par- 
ticle through. The escort "informs" the particle 
about the environment. Such an insight consid- 
ering behavior of the particle with the point of 
view of Feynman path integrals, can be more 
productive, than fantastic ideas about splitting 
particles passing through the slit grating and 
their confluence as soon as the slits are left far 
behind. 

We need in this connection to continue discus- 
sion of the problems, relating to the virtual and 
real trajectories. Let us consider a number of 
variational procedures applied to the complex- 
valued function ip{q,p;t), wave function. It per- 



mits to see which variational scheme ends up 
the classical equations of movement and what 
scheme leads to the quantum mechanical equa- 
tions. The problem, in fact, is to retrace how the 
virtual trajectories relate to the real trajectories, 
Bohmian trajectories. 

A. Classical domain 

Let the complex-valued function be as follows 



(34) 



= v^- exp<j ^ / L(g,g;r)rfr 



to 



Let us demand that this function would retain a 
constant value along a path from the initial time 
t = to to final t = ti. In other words, varia- 
tion of this function along this time interval has 
to vanish, i.e., Sip = 0. Applying this varia- 



tion to the expression (34), we disclose, that the 



complex-valued function ends by two equations, 
separately for real and imaginary parts. Each 
equation should vanish. 

Real part of the above equation leads imme- 
diately to the continuity equation ^ for the 
probability density p. Imaginary part reduces 
to variation of the action S, see equation ([T]), 
along paths from the initial time to to the final 
time ti- After a series of mathematical transfor- 
mations [1] we disclose, that the particle moves 
along an optimal path that is described by the 
Hamilton- Jacobi equation 

We observe in the classical case, that the parti- 
cles move along the classical trajectories submit- 
ting to the principle of least action. By moving 
along the optimal paths, ensemble of the par- 
ticles resembles a "cloud" having a density dis- 
tribution p. Evolution of the cloud obeys the 
continuity equation for the density distribution. 

B. Quantum domain 

In contrast to the previous searching of a single 
trajectory connecting the initial and final points, 
here all trajectories connecting these points are 
to be considered. They pass through all inter- 
mediate points belonging to a conditional set 
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TZ^. And these all paths have to be evaluated 
jointly. Such a description goes back to the in- 
teger Chapman-Kolmogorov equation 



z;t + s) = J p{x, y; t)p{y, z; s)dy (35) 

which gives transitional probability densities of 
a Markov sequence. 




FIG. 16: Computation of all possible paths passing from 
point go to point q\ through possible intermediate points q^c £ 
TZ^ represents a core of the path integral method. 



Essential difference from the classical proba- 
bility theory is that instead of the probabilities 
quantum mechanics deals with amplitudes con- 
taining imaginary terms, that bear phase infor- 
mation. The amplitude function charged with a 
phase term is a complex-valued function called 
the wave function. In this key, a transition from 
an initial state go to a final state qi through all 
intermediate instances qx given on a conditional 
set 7^^, see Fig. 9, is represented by the following 
path integral 

?/'(gi,go,i + 5t) = 
K{qi, q^; t + 6t, t)ip{qx, go; t)V^qx (36) 



7^3 

in the limits St ^ and qi ^ qx- Here symbol 
T>^qx represents a differential element of volume 
in the set TZ^. Integral kernel, i.e., propagator, 
K{qi,qx;t+6t,t) we suppose has a standard form 

K{qi,qx;t+St,t) = ^ exp|^L(g;, g;)5t|. (37) 

Standard Lagrangian in this place has a form [7] 



^(g-;,i) = fP^ -Uiqx). (38) 



Here U{qx) is a potential energy of the particle 
localized at the point qx G 71^. And (g! — qx)/St 
is a velocity g^ attached to the same point g^ 
and oriented in the direction of the point gl. 
The next step is to expand terms, ingoing into 



the integral (36), into Taylor series. For example 



the wave function written at the left is expanded 
up to the first term 

HQu go, t + 5t)^ ^(gl, go, t) + -^5t. (39) 

As for the terms under the integral, here we pre- 
liminarily make some transformations. Let us 
determine a small increment 



(40) 



It is believed, that a main contribution in the 
integral is only given by intermediate sources 
from the set TZ^ situated near the point gl, i.e., 
corresponding to small ^. In this case, the La- 



grangian (38) can be rewritten as 



Liqx,qx) 



2 6t 



2-Wi-O, (41) 



where the the potential energy f/(gi — is sub- 
jected to expansion into the Taylor series by the 

— * 

small parameter ^. The under integral wave 
function ipiQx, qo', t) = i'iqi — ^, go; t) is subjected 
also to expansion into the Taylor series up to the 
second terms of the expansion 



(42) 



Taking into account the expressions (39 )-42 ) and 



substituting theirs into Eq. 36 ) we get 

i m 



, , ^ ^ \ dip ^ 
V^(gi,go,t) + -^dt 
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n 

gl, g-o, t) - (V^ ■ O + ■ e/'2)v% (43) 



Here the exponent, containing the potential term 
U (gi —^)St, was preliminarily expanded into the 
Taylor series. The term in this expansion having 
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a factor of the third order of smallness, in 
what follows is ignored. 

The free term ip{qi,qo,t) represented both 
from the left and from the right can disappear, 
if and only if an expression, containing the con- 
stant A, will satisfy the following condition 

(44) 

from here it is follows 



A = - 



/27riMt\ 



3/2 



The power 3 emerges here because that the inte- 
gration is fulfilled on the 3-dimensional set 'R? . 
In addition to this observation, we need to clar- 
ify integration of the terms (VV' • ^ ) and A?/) ■ ^ . 
With this aim in the mind, we mention the fol- 
lowing two integrals [7] 

1 /■ {\m^^ 



A 



exp-! 



TV- 



h 2 6t 



>^V^^ = 0. 



and 



1 

A 



exp 



7^3 




-6t. 



m 



(46) 



(47) 



In accordance with the first integral, contribu- 
tions of the terms V0 and Vf/ in the expres- 



sion (43) disappear. Whereas, the second term. 



A'?/', gains the factor (\h6t/m)/2. 
Taking into account the above stated, we re- 



turn to Eq. (43) and write out only terms, that 



have the factors of the first order of smallness 
regarding 8t 



ot 



■—Uiqi)iljiqi,qo,t) + —A^p. 



(4J 



So far, we retained here the free term ip{qi, %, t). 
Meanwhile, other terms, having orders of small- 
ness 5t^ and smaller, are excluded from this 
equation. Finally, by reducing ip{qi,qo,t) from 
the left and from the right, we come to the 
Schrodinger equation 



ih 



dip 



2m 



Aip + U{q)ip, 



(49) 



describing the wave field ip in the configuration 
set TZ^. This set evolves by force of infinitesimal 
transformations {C, — )■ 0) by shifting in the time 
{6t 0), see Fig. 
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Obviously, the subscript 1 
now can be dropped. 

It should be noted, since at derivation of 
the Schrodinger equation we use a condition of 
smallness \qi — qx\ ^ 1, then decay of radiation 
by secondary sources with distance is not take 
into account. That is, in this derivation only 
a small set of the secondary sources, nearest to 
the point qi, give real contribution to the path 



(45) integral (36). The distant sources give so small 



contribution, that we neglect them. 
Now, let us substitute into Eq. (49) the wave 



function ip to be represented by multiplication of 
the amplitude part p^/^ on the phase exponent 
exp{iS/h}, see Eq. (34). Next, by separating 



the imaginary and real parts of the Schrodinger 
equation [28] we come to equations (27)-(28). 



They are the quantum HJ-equation and the con- 
tinuity equation, respectively. 



VII. 



CONCLUSION: VIRTUAL AND REAL 
TRAJECTORIES 



Let us return to Fig. 16 Spherical waves out- 



lined in this Figure demonstrate radiation from 
each point q^ through which virtual trajectories 
can pass. Their wavelengths are congruous with 
that of the particle. Those spherical waves radi- 
ated from each point G 7^3, reaching gl, create 
in this point an interference effect. This effect 
shows whether a real particle passes through the 
point qi. If this event can happen, though, what 
kind is chance of it? 

One can see, that we outlined here the 
Huygens-Fresnel principle [22] ■ It proclaims that 
each point q^ at an advanced wave front is in fact 
the center of a fresh disturbance and the source 
of a new wave radiation. And the advancing 
wave as a whole may be regarded as the sum 
of all the secondary waves arising from points in 
the medium already traversed. All the secondary 
waves are coherent, since they are activated from 
the one source given in go- 

It is important to emphasize, that all rays from 
such secondary sources represent in fact virtual 
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trajectories passing from the sources up to the 
point qi. Along with the other virtual trajecto- 
ries generated by the other secondary sources, all 
together they create in the point qi an averaged 
effect of contribution of these secondary sources. 
This averaged effect shows whether a real parti- 
cle passes by this route and what probability of 
this event can be. 

One can imagine, that all physical space is 
filled by such secondary sources of spherical 
waves. These sources are virtual in that sense, 
that they manifest themselves via zero-point os- 
cillations of vacuum. It means, that there are 
permanent creations and annihilations of virtual 
pairs "particle-antiparticle". Lifetime of such 
pairs are associated with their total energy by 
the Heisenberg uncertainty principle, i.e., the 
larger uncertainty of the total energy the smaller 
the lifetime of the pair. The same can be said 
about their scattering from each other: if the po- 
sition is known perfectly, then the momentum is 
completely unknown, and vice versa. 

Imagine now, let an area of scattering of such 
pairs be equal to about an effective width of the 
slit, 0", that is multiple of the wavelength. In 
that case, the momentum, p = h/a, can be suffi- 
cient to force the particle to change a direction of 
scattering. Virtual spherical waves in the vicin- 
ity of the slit can effect to behavior of the particle 



passing through the slit. Ensemble of such vir- 
tual waves, each attributed to its slit, creates a 
polarized effect similar to the domino effect. In 
other words, in a space, adjoining to the slit grat- 
ing, zero-point oscillations of vacuum are quite 
ordered. The slits in effect control the vacuum 
fluctuations. It is not analogous to ordered po- 
sitions of atoms in crystals. The ordering of vir- 
tual oscillations of vacuum, in the case of the 
slit grating, is expressed via organizing the wave 
function that describes transport of the particles 
through the grating. The wave function repre- 
sents itself the de Broglie pilot wave [SUl EI] that 
guides a particle from a source to a detector. A 
path from the source to the detector, called the 
Bohmian trajectory, is resulted from solution of 
the guidance equation [T9] . 
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